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Abstract 

Recent results on particle momentum and spin correlations are discussed in view 
of the role played by the effects of quantum statistics, including multiboson and 
coherence phenomena, and final state interaction. Particularly, it is demonstrated 
that the latter allows for (i) correlation femtoscopy with unlike particles; (ii) study 
of the relative space-time asymmetries in the production of different particle species 
(e.g., relative time delays or spatial shifts due to collective flows); (iii) study of the 
particle strong interaction hardly accessible by other means (e.g., in AA system) 
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1 Introduction 



The momentum correlations of particles at small relative velocities are widely used to 
; study space-time characteristics of the production processes, so serving as a correlation 
femtoscope. Particularly, for non-interacting identical particles, like photons or, to some 
extent, pions, these correlations result from the interference of the production amplitudes 
due to the symmetrization requirement of quantum statistics (QS) [1, 2}. There exists 
[3] a deep analogy of the momentum QS correlations of photons with the space-time 
correlations of the intensities of classical electromagnetic fields used in astronomy to 
measure the angular radii of stellar objects based on the superposition principle - so 
"q ■ called HBT intensity interferometry [4]. 2 

The momentum QS correlations were first observed as an enhanced production of the 
pairs of identical pions with small opening angles (GGLP effect [1]). Later on, Kopylov and 
Podgoretsky [2] settled the basics of correlation femtoscopy; particularly, they suggested 
to study the interference effect in terms of the correlation function and clarified the role 
of the space-time characteristics of particle production in various physical situations. 

The momentum correlations of particles emitted at nuclear distances are also influ- 
enced by the effect of final state interaction (FSI) [6, 7, 8]. Thus the effect of the Coulomb 
interaction dominates the correlations of charged particles at very small relative momenta 
(of the order of the inverse Bohr radius of the two-particle system), respectively suppress- 
ing or enhancing the production of particles with like or unlike charges. Though the FSI 

1 Work supported by grant 202/01/0779 of the Grant Agency of the Czech Republic. 
2 This analogy is sometimes misunderstood and the momentum correlations are mixed up with the 
space-time HBT correlations in spite of their orthogonal character and the failure of the superposition 
principle for correlations of identical fermions. In fact, in spite of the common QS origin of the momentum 
correlations of identical particles and the space-time HBT correlations (allowing for a generalization of the 
latter to any type of identical bosons or fermions), the corresponding correlation measurements differ in 
principle [3] (see also [5]). The former, being the momentum-energy measurement, yields the space-time 
picture of the source, while the latter does the opposite. In particular, the dependence of the number 
of coincident two-photon counts on the distance between detectors (a quantum analogy of the HBT 
measurement) provides the information on the characteristic relative three-momenta of emitted photons 
and so, when divided by the mean detected momentum, on the angular size of a star but, of course, - no 
information on the star radius or lifetime. 
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effect complicates the correlation analysis, it is an important source of information allow- 
ing for the coalescence femtoscopy (see, e.g., [9, 10, 11, 12]), the correlation femtoscopy 
with unlike particles [8, 13] including the access to the relative space-time asymmetries 
in particle production [14] and a study of particle interaction hardly accessible by other 
means. 

We do not touch here the fluctuation measures which are closely related with particle 
correlations in momentum space and carry an important information on the dynamics 
and space-time evolution of the production process (see [15] for a recent review). 

The rest of the report is organized as follows. In section 2, we briefly review the formal- 
ism of particle correlations at small relative velocities. The basic concepts of femtoscopy 
with identical and nonidentical particles, including the access to the relative space-time 
shifts in the emission of various particle species, and some recent results are reviewed in 
sections 3, 5 and 7. In section 4, we discuss the present theoretical and experimental 
status of the multiboson and coherence phenomena in multiparticle production. Recent 
results from correlation measurements of the strong interaction in various two-particle 
systems are reviewed in section 6. In section 8, we briefly discuss spin correlations as a 
new femtoscopy tool. We conclude in Section 9. 

2 Formalism 

The ideal two-particle correlation function TZ(p\ ) p2) is defined as a ratio of the measured 
two-particle distribution to the reference one which would be observed in the absence of 
the effects of QS and FSI. In practice, the reference distribution is usually constructed 
by mixing the particles from different events of a given class, normalizing the correlation 
function to unity at sufficiently large relative velocities. 

Usually, it is assumed that the correlation of two particles emitted with a small relative 
velocity is influenced by the effects of their mutual QS and FSI only 3 and that the mo- 
mentum dependence of the one-particle emission probabilities is inessential when varying 
the particle four-momenta pi and P2 by the amount characteristic for the correlation due 
to QS and FSI (smoothness assumption) . Clearly, the latter assumption, requiring the 
components of the mean space-time distance between particle emitters much larger than 
those of the space-time extent of the emitters, is well justified for heavy ion collisions. 

The correlation function is then given by a square of the properly symmetrized Bethe- 
Salpeter amplitude in the continuous spectrum of the two-particle states, averaged over 
the four-coordinates Xi = {tj, r^} of the emitters and over the total spin S of the two- 
particle system [8]. After the separation of the unimportant phase factor due to the c.m.s. 
motion, this amplitude reduces to the one depending only on the relative four-coordinate 
Ax = x\ — X2 = {t,r} and the generalized relative momentum q = q — P(qP)/P 2 , where 
q = pi — p 2 , P = pi + P2 and qP = mi 2 — m 2 2 ; in the two-particle c.m.s., P = 0, 
q = {0,2k*} and Ax = {t*,r*}. At equal emission times of the two particles in their 
c.m.s. (t* = t* — = 0), the reduced non-symmetrized amplitude coincides with a 
stationary solution (r*) of the scattering problem having at large distances r* the 
asymptotic form of a superposition of the plane and outgoing spherical waves (the minus 

3 Besides the events with a large phase-space density fluctuations, this assumption may not be justified 
also in low energy heavy ion reactions when the particles are produced in a strong Coulomb field of residual 
nuclei. To deal with this field a quantum adiabatic (factorisation) approach can be used [16]. 



2 



sign of the vector k* corresponds to the reverse in time direction of the emission process). 
The Bethe-Salpeter amplitude can be usually substituted by this solution (equal time 
approximation) 4 so, for non-identical particles, 

n( Pl , P2 ) = Y,Ps(\&(r*)\ 2 )s; (i) 

s 

for identical particles, the amplitude in Eq. (1) enters in a symmetrized form: 

VdtV) - hffiV) + (-l) 5 l£ (+) (r')]/V2. (2) 

The averaging in Eq. (1) is done over the four-coordinates of the emitters at a given total 
spin S of the two-particles, ps is the corresponding population probability, Y,s Ps — 1- F° r 
unpolarized particles with spins s 1 and s 2 the probability ps = (2S , + l)/[(2s 1 + l)(2s 2 + l)]- 
Generally, the correlation function is sensitive to particle polarization. For example, if 
two spin-1/2 particles are initially emitted with polarizations V\ and 7*2 then [8] 

p = (l-7VP 2 )/4, p 1 = (3 + 7VP 2 )/4. (3) 



3 Femtoscopy with identical particles 

For identical pions or kaons, the effect of the strong FSI is usually small and the effect of 
the Coulomb FSI can be in first approximation simply corrected for (see [17] and references 
therein). The corrected correlation function is determined by the QS symmetrization only 
(see Eq. (2) and substitute the non-symmetrized amplitude by the plane wave e iqx ^ 2 ): 

ll(p 1 ,p 2 ) = l + (cos(qAx)). (4) 

Its characteristic feature is the presence of the interference maximum at small components 
of the relative four-momentum q with the width reflecting the inverse space-time extent of 
the effective production region. For example, assuming that for a fraction A of the pairs, 
the pions are emitted independently according to one-particle amplitudes of a Gaussian 
form characterized by the space-time dispersions r 2 and Tq while, for the remaining frac- 
tion (1 — A) related to very long-lived sources (77, r/, K®, A, . . . ), the relative distances 
r* between the emitters in the pair c.m.s. are extremely large, one has 

KipuPi) = 1 + A exp (-r 2 q 2 - r 2 g 2 ) = 1 + A exp (-r 2 q^ - (r 2 + v 2 r 2 )q 2 L ) , (5) 

where qx and qi are the transverse and longitudinal components of the three-momentum 
difference q with respect to the direction of the pair velocity v = P/Po- One may see that, 
due to the on-shell constraint [2] q = vq = vq^ (following from the equality qP = 0), 
strongly correlating the energy difference go with the longitudinal momentum difference 

4 For non-interacting particles, the non-symmetrized Bethe-Salpeter amplitude reduces to the plane 
wave e l Q x / 2 = e~* k * r * which is independent of the relative time in the two-particle c.m.s. and so, coincides 
with the corresponding equal-time amplitude. For interacting particles, the equal time approximation is 
valid on condition [8] \t*\ <C m 2j ir* 2 for sign(t*) = ±1 respectively. This condition is usually satisfied 
for heavy particles like kaons or nuclcons. But even for pions, the t* — approximation merely leads 
to a slight overestimation (typically < 5%) of the strong FSI effect and, it doesn't influence the leading 
zero-distance (r* <§C \a\) effect of the Coulomb FSI. 
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qL, the correlation function at vt > r substantially depends on the direction of the 
vector q even in the case of a spherically symmetric spatial form of the production region. 

Note that the on-shell constraint makes the g-dependence of the correlation function 
essentially three-dimensional (particularly, in pair c.m.s., qAx = — 2k*r*) and thus makes 
impossible the unique Fourier reconstruction of the space-time characteristics of the emis- 
sion process. However, within realistic models, the directional and velocity dependence 
of the correlation function can be used to determine both the duration of the emission 
and the form of the emission region [2], as well as - to reveal the details of the production 
dynamics (such as collective flows; see, e.g., [18, 19] and the reviews [20, 21]). For this, 
the correlation functions can be analyzed in terms of the out (x), side (y) and longitudinal 
(z) components of the relative momentum vector q = {q x ,q y ,qz} [22, 23]; the out and 
side denote the transverse, with respect to the reaction axis, components of the vector 
q, the out direction is parallel to the transverse component of the pair three-momentum. 
The corresponding correlation widths are usually parameterized in terms of the Gaussian 
correlation radii Ri, 

K(p 1 , P2 ) = 1 + A eM-Rlll - R 2 y Q 2 y - RW Z ~ RlzQxQz) (6) 

and their dependence on pair rapidity and transverse momentum is studied. The form of 
Eq. (6) assumes azimuthal symmetry of the production process [20, 22]. Generally, e.g., 
in case of the correlation analysis with respect to the reaction plane, all three cross terms 
qiqj contribute [24]. 

It is well known that particle correlations at high energies usually measure only a 
small part of the space-time emission volume, being only slightly sensitive to its increase 
related to the fast longitudinal motion of particle sources. In fact, due to limited source 
decay momenta p^ of few hundred MeV/c, the correlated particles with nearby velocities 
are emitted by almost comoving sources and so - at nearby space-time points. In other 
words, the maximal contribution of the relative motion to the correlation radii in the two- 
particle c.m.s. is limited by the moderate source decay length rp^/m. The dynamical 
examples are sources-resonances, colour strings or hydrodynamic expansion. To substan- 
tially eliminate the effect of the longitudinal motion, the correlations can be analyzed in 
terms of the invariant variable qw, = Q — (—q 2 ) 1 ^ 2 = 2k* and the components of the 
momentum difference in pair c.m.s. (q* = Q = 2k*) or in the longitudinally comoving 
system (LCMS) [25]. In LCMS each pair is emitted transverse to the reaction axis so 
that the generalized relative momentum q coincides with q* except for the component 
q x = ltq* x i where j t is the LCMS Lorentz factor of the pair. 

Particularly, in the case of one-dimensional boost invariant expansion, the longitudinal 
correlation radius in the LCMS reads [19] R z ~ (T/m t ) 1 / 2 r, where T is the freeze-out 
temperature, r is the proper freeze-out time and m t is the transverse particle mass. In 
this model, the side radius measures the transverse radius of the system while, similar 
to Eq. (5), the square of the out radius gets an additional contribution (p t /mt) 2 Ar 2 
due to the finite emission duration At. The additional transverse expansion leads to a 
slight modification of the p t -dependence of the longitudinal radius and - to a noticeable 
decrease of the side radius and the spatial part of the out radius with p t . Since the freeze- 
out temperature and the transverse flow determine also the shapes of the m t -spectra, the 
simultaneous analysis of correlations and single particle spectra for various particle species 
allows to disentangle all the freeze-out characteristics (see the review [20]). It appears 
that with the increasing energy of heavy ion collisions from AGS and SPS up to the 
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highest energies at RHIC, the data show rather weak energy dependence [26] and point 
to the kinetic freeze-out temperature somewhat below the pion mass, a strong transverse 
flow (with the mean transverse flow velocity at RHIC exceeding half the velocity of light 
[27]), a short evolution time of 8-10 fm/c and a very short emission duration of about 
2 fm/c. The short evolution and emission duration at RHIC are also supported by the 
correlation analysis with respect to the reaction plane [28]. The small time scales at RHIC 
were not expected in transport and hydrodynamic models [29, 30] and may indicate an 
explosive character of particle production (see, e.g., [31, 32]). In fact, the RHIC data can 
be described in so called blast wave model [33, 34] assuming a strong three-dimensional 
expansion with a sharp boundary of the freeze-out density profile in transverse plane. The 
same model with ~ 15% lower mean transverse flow velocity is also consistent with the 
SPS data [35]. 

4 Multiboson and coherence effects 

In present and future heavy ion experiments at SPS, RHIC and LHC many hundreds 
or thousands of pions can be produced per a unit rapidity interval. Since pions are 
bosons there can be multiboson effects enhancing the production of pions with low relative 
momenta thus increasing the pion multiplicities, softening their spectra and modifying the 
correlation functions (see [36, 37, 38] and references therein). In particular, it was shown 
[37] that the width of the low-pt enhancement due to BE condensation decreases with 

— 1/2 

the system size as r and this narrowing makes easier the identification of this effect 
among others. For the events of approximately fixed multiplicity, the multiboson effects 
can be triggered by decreasing correlation strength and a dip in the two-pion correlation 
function at intermediate relative momenta [37, 38]. 

Though the present data does not point to any spectacular multiboson effects, one 
can hope to observe new interesting phenomena like boson condensation or speckles in 
some rare events or in eventually overpopulated kinematic regions with the pion density in 
the 6-dimensional phase space, / = (27r) 3 d 6 n/<i 3 p<i 3 x, of the order of unity. An example 
is a rapidly expanding system with the entropy much smaller than in the case of total 
equilibrium. Then a strong transverse flow can lead to rather dense gas of soft pions in 
the central part of the hydrodynamic tube at the final expansion stage (see, e.g., [39]). 
Another reason can be the expected formation of quark-gluon plasma or mixed phase. Due 
to large gradients of temperature or velocity the hydrodynamic layer near the boundary 
with vacuum can decay at a large phase space density and lead to pion speckles even at 
moderate transverse momenta [40]. 

In the low-density limit (/ <C 1), the mean phase space density at a given momentum p 
can be estimated as the mean number of pions interfering with a pion of momentum p (ra- 
pidity y and transverse momentum p t ) and building the Bose-Einstein (BE) enhancement 
in the two-pion correlation function [41, 42]: (/) p ~ 7T 3 / 2 N(p)/V, where N(p) = d 3 n/d 3 p 
and V = r x r y r z is the interference volume defined in terms of the outward (r x ), sideward 
(r y ) and longitudinal (r z ) interferometry radii. Typically (/) p ~ 0.1 for mid-rapidities 
and p t ~ (pt) [41]. The data are also consistent with the phase space density of pions 
near the local thermal equilibrium [43, 44]. 

At AGS and SPS energies the interference volume V seems to scale with dn/dy (see, 
e.g., [44, 45]) pointing to the freeze-out of pions at a constant phase space density. This 
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trend is however questioned by recent STAR data from RHIC, indicating an increase of 
the freeze-out phase space density with energy (a slight increase of V is not sufficient to 
balance ~ 50% increase of dn/dy as compared with SPS) and centrality [46]. Extrapola- 
tion of the RHIC phase space density measurements to low transverse momenta predicts 
(/) p close to unity for central events, suggesting that significant multiboson effects can 
be present at low p t at RHIC. 

According to lattice Monte Carlo calculations including dynamical fermions, decon- 
fining phase transition leading to a quark-gluon plasma (QGP) phase of matter is accom- 
panied by restoration of chiral symmetry. Subsequent phase transition into the hadronic 
phase can be revealed, particularly, through substantial delays in particle emission and/or, 
through the coherent component of the pion radiation. This component would be charac- 
terized by a narrow Poisson multiplicity distribution, contrary to wide multiplicity fluc- 
tuations in the usual BE condensate. The pions in the coherent state may appear from 
the decay of a quasi-classical pion field (the order parameter of the phase transition), the 
latter possibly related to the spontaneous chiral symmetry breaking via the formation of 
the disoriented chiral condensate (DCC) (see [47] and a review [48]). 

The most plausible mechanism of DCC formation is a fast expansion of hot QGP 
resulting in a rapid supression of thermal fluctuations (quenching), which in turn triggers 
a dramatic amplification of soft pion modes. The detection and study of DCC is expected 
to provide valuable information about the chiral phase transition and vacuum structure 
of strong interactions. DCC formation is usually expected to be associated with large 
event-by-event fluctuations in the ratio of neutral to charged pions in a certain phase- 
space domain. The search for these fluctuations at CERN SPS has so far resulted in 
setting only an upper limit on the production of a single DCC domain [49]. The absence 
of experimental evidence for isospin fluctuations has been however recently claimed to be 
in agreement with presumably more realistic picture of an " unpolarized" DCC with the 
Fourier modes of the field randomly oriented in isospin space (instead of being aligned 
as in the original DCC scheme) [50]. The search for other DCC signatures like low 
momentum pion clusters is therefore important. Particularly, one can exploit the impact 
of the admixture of coherent radiation on the QS and Coulomb correlations of like and 
unlike pions [51]. Other possibilities of experimental investigations of BE condensate and 
DCC phenomena have been discussed, e.g., in [52, 53]. 

The presence of the coherent pions (or pions emitted in the same quantum state) 
manifests itself also as a suppression of the BE correlations of two or more identical pions 
[7, 54, 55, 56]. Unfortunately, there are also other reasons leading to the suppression of 
particle correlations. Besides the experimental effects like finite resolution and particle 
misidentification (that can be corrected for), presumably the most important one is the 
contribution of the particles emitted by long-lived sources [57], leading to the appearance 
of the parameter A < 1 in Eqs. (5) and (6). Also the usual Gaussian parameterizations 
of the QS correlation functions may be inadequate and lead to A < 1 in the presence of 
the sources with moderate but very different space-time characteristics [57, 58, 59]. 

In principle, the effect of long-lived sources can be eliminated in a combined analysis 
of two-pion and three-pion correlation functions. The measured quantity is the genuine 
three-pion correlation normalized with the help of the three two-pion contributions - its 
intercept measures the chaotic or coherent fraction [60]. First such measurements have 
been done only recently in heavy ion experiments at CERN SPS [61, 62] and RHIC [63] 
and, in e + e~ collisions at LEP [64]. The most accurate ones at RHIC and LEP indicate 
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a dominant chaotic fraction though the systematic errors allow for a substantial coher- 
ent component. Some sources of the systematic errors, e.g., the simplified treatment 
of the two-body Coulomb and strong FSI, can be overcome. However others, e.g., the 
approximate (factorization) treatment of the multiparticle FSI or the insufficiently dif- 
ferential analysis of the three-pion correlation function, can hardly be avoided at present 
computational and experimental possibilities. 

5 Femtoscopy with unlike particles 

The complicated dynamics of particle production, including resonance decays and particle 
rescatterings, leads to essentially non-Gaussian tail of the distribution of the relative 
distances r* of the particle emitters in the pair rest frame. Therefore, due to different 
r*-sensitivity of the QS, strong and Coulomb FSI effects, one has to be careful when 
analyzing the correlation functions in terms of simple models. Thus, the QS and strong 
FSI effects are influenced by the r*-tail mainly through the suppression parameter A 
already for distances of the order of inverse g-resolution (typically some tens fm) while, 
the Coulomb FSI is sensitive to the distances as large as the pair Bohr radius |a|; for 
7T7T, irK, Tip, KK, Kp and pp pairs, \a\ = 387.5, 248.6, 222.5, 109.6, 83.6 and 57.6 fm, 
respectively. Clearly, the usual Gaussian parameterizations of the distributions of the 
components of the distance vector r* may lead to inconsistencies in the treatment of QS 
and FSI effects (the Coulomb FSI contribution requiring larger effective radii). These 
problems can be at least partially overcome with the help of transport code simulations 
accounting for the dynamical evolution of the emission process and providing the phase 
space information required to calculate the QS and FSI effects on the correlation function. 

Thus, in a preliminary analysis of the NA49 correlation data from central Pb + Pb 158 
AGeV collisions [65, 66], the freeze-out phase space distribution has been simulated with 
the RQMD v. 2. 3 code [67]. The correlation functions have been calculated using the code 
of Ref. [8], weighting the simulated pairs by squares of the corresponding wave functions. 
The dependence of the correlation function on the invariant relative momentum Q = 2k* 
was than fitted according to the formula [65] 

7Z(Q) = norm • [purity • RQMD(r* — > scale • r*) + (1 — purity)]; (7) 

to account for a possible mismatch in (r*), the dependence on the r*-scale parameter has 
been introduced using the quadratic interpolation of the points simulated at three scales 
chosen at 0.7, 0.8 and 1. The fitted values of the purity parameter are in reasonable 
agreement with the expected contamination of ~ 15% from strange particle decays and 
particle misidentification. The fitted values of the scale parameter indicate that RQMD 
overestimates the distances r* by 10-20%. Similar overestimation has been also observed 
when comparing RQMD predictions with the NA49 data on pp and 7r ± 7r ± correlations 
[68, 69, 70]. 

Recently, there appeared data on pA correlation functions from Au + Au experiment 
E985 at AGS [71] and Pb + Pb experiment NA49 at SPS CERN [72]. As the Coulomb FSI 
is absent in pA system, one avoids here the problem of its sensitivity to the r*-tail. Also, 
the absence of the Coulomb suppression of small relative momenta makes this system 
more sensitive to the radius parameters as compared with pp correlations [73] . In spite of 
rather large statistical errors, a significant enhancement is seen at low relative momentum, 
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consistent with the known singlet and triplet pA s-wave scattering lengths. In fact, the 
fits using the analytical expression for the correlation function (originally derived for pn 
system [8]) yield for the AGS data [66] the purity of 0.5 ± 0.2 and the Gaussian radius of 
4.5 ± 0.7 fm. For the NA49 data the fitted parameters are [72] 0.17 ± 0.11 and 2.9 ± 0.7 
fm. The fitted AGS purity is consistent with the estimated one, while the NA49 purity 
is about one standard deviation too low. Fixing the NA49 purity at the estimated value 
of 0.33, the Gaussian radius increases by about 1 fm and becomes 3.8 ± 0.4 fm [72]. The 
fitted AGS and NA49 radii are in agreement with the radii of 3-4 fm obtained from pp 
correlations in heavy ion collisions at GSI, AGS and SPS energies. 

6 Correlation measurement of strong interaction 

In case of a poor knowledge of the two-particle strong interaction, which is the case for 
meson-meson, meson-hyperon or hyperon-hyperon systems, 5 it can be improved with the 
help of correlation measurements. 

In heavy ion collisions, the effective radius tq of the emission region can be considered 
much larger than the range of the strong interaction potential. The FSI contribution is 
then independent of the actual potential form [75]. At small Q = 2k*, it is determined 
by the s-wave scattering amplitudes f s (k*) [8]. In case of \f s \ > r , this contribution is 
of the order of \f s /r \ 2 and dominates over the effect of QS. In the opposite case, the 
sensitivity of the correlation function to the scattering amplitude is determined by the 
linear term f s /r . 

The possibility of the correlation measurement of the scattering amplitudes has been 
demonstrated [66] in a recent analysis of the NA49 it + it~ correlation data within the 
RQMD model. For this, the strong interaction scale has been introduced (similar to the r*- 
scale), redefining the original s-wave 7r + 7r~scattering length f = 0.232 fm: f Q — > sisca- /o- 
The fitted parameter sisca = 0.63 ± 0.08 appears to be significantly lower than unity. To 
a similar shift (~ 20%) point also the recent BNL data on decays [76]. These results 
are in agreement with the two-loop calculation in the chiral perturbation theory with a 
standard value of the quark condensate [77]. 

Recently, also the singlet A A s-wave scattering length fo has been estimated [66, 72] 
based on the fits of the NA49 AA data. Using the analytical expression for the correlation 
function [78] (originally derived for nn system [8]) and fixing the purity of direct A- 
pairs at the estimated value of 0.16 and varying the effective radius r in the acceptable 
range of several fm, one gets [72] e.g., f = 2.4 ± 2.1 and 3.2 ± 5.7 fm for r = 2 and 
4 fm respectively (we use the same sign convention as for meson-meson and meson- 
baryon systems). Though the fit results are not very restrictive, they likely exclude the 
possibility of a large positive singlet scattering length comparable to that of ~20 fm for 
the two-nucleon system. 

The important information comes also from AA correlations at LEP [79]. Here the ef- 
fective radius r is substantially smaller than the range of the strong interaction potential, 
so the AA correlation function is sensitive to the potential form and requires the account 
of the waves with orbital angular momentum up to Z ~ 20 [80]. In Ref. [79], the strong 
interaction has been neglected and the observed decrease of the AA correlation function 

5 The AA system is of particular interest in view of an experimental indication on the enhanced AA 
production near threshold [74] and its possible connection with the 6-quark H dibaryon problem. 
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at small Q has been attributed solely to the effect of the QS (Fermi-Dirac) suppression. 
The correlation function has been fitted by the expression 6 

-R=l- l -\{l + V 2 )eM-rlQ 2 ) (8) 

corresponding to the simple Gaussian distribution of the components of the relative dis- 
tance vector r* characterized by a dispersion 2r%. The fit results are however unsatisfac- 
tory for two reasons [80]: (i) the parameter A = 1.2 ± 0.2 (neglecting in Eq. (8) the "P 2 
polarization term on a percent level) is significantly higher than the value of ~ 0.5 ex- 
pected due to the feed-down from E° and weak decays; (ii) the parameter r = 0.11 ±0.02 
fm appears to be smaller than the string model lower limit of ~ 0.2 fm. Therefore, the 
observed anti-correlation at small Q can be considered as a direct evidence for a repul- 
sive core in the AA interaction potential. 7 In fact, reasonable fits can be achieved using 
the Nijmegen singlet potential NSC97e [81], rescaling the triplet one from Ref. [82] and, 
neglecting spin-orbit and tensor couplings. For example, at a fixed A = 0.6, the fitted 
radius takes an acceptable value r = 0.29 ± 0.03 fm [80]. 

7 Accessing relative space— time asymmetries 

The correlation function of two non-identical particles, compared with the identical ones, 
contains a principally new piece of information on the relative space-time asymmetries 
in particle emission such as mean relative time delays in the emission of various particle 
species [14]. It can be particularly useful in searches for the effects of the quark-gluon 
plasma phase transition like delays between the emission of strange and antistrange par- 
ticles due to the process of strangeness distillation from the mixed phase. The important 
information is contained also in the spatial part of the asymmetry related, in particular, 
with the intensity of the collective flow [66]. 

Since the information on the relative space-time shifts enters in the two-particle wave 
function through the terms odd in k*r* = p^(r* — r^), it can be accessed studying the 
correlation functions lZ +i and 1Z_i with respectively positive and negative projection k* 
of the momentum k* = = — p*. on a given direction i or, - the ratio lZ + i/lZ-i. For 
example, i can be the direction of the pair velocity or, any of the out (x), side (y), 
longitudinal (z) directions. Note that in the LCMS system, 

r* = Ax* = lt (Ax - v t At), r* y = Ay* = Ay, r* z = Az* = Az, (9) 

where ^ t = (1 — v t 2 ) 1 / 2 and v t = P t /Po are the pair LCMS Lorentz factor and velocity. 
One may see that the asymmetry in the out (x) direction depends on both space and 
time asymmetries (Ax) and (At). In case of a dominant Coulomb FSI, the intercept of 
the correlation function ratio is directly related with the asymmetry (r*) [83, 84] (see also 
[85]): 

ft+i/ft-i » 1 + 2(r*)/a, (10) 

6 The singlet and triplet contributions to the correlation function 1Z = 1Z S + IZt are 1Z s ,t — Ps,t[^ ^ 
A exp(— TqQ 2 )], where p St t depend on the A-polarization T according to Eq. (3) with "Pi = Vi = "P- 

7 The repulsive core arises due to the exchange of vector mesons and is present, e.g., in various Nijmegen 
potentials used for the analysis of the double A hypernuclei. The core height and width are about 9 GeV 
and 0.4 fm respectively. The s-wave scattering length (effective radius) ranges from about 0.3 (15) fm to 
11 (2) fm. 
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where a = (fiziz 2 e 2 )~ 1 is the Bohr radius of the two-particle system taking into account 
the sign of the interaction (^e are the particle electric charges, /i is their reduced mass). 

At low energies, the particles in heavy ion collisions are emitted with the characteristic 
emission times of tens to hundreds fm/c so that the observable time shifts should be of 
the same order [14]. Such shifts have been indeed observed with the help of the 1Z+/1Z- 
correlation ratios for proton-deuteron systems in several heavy ion experiments at GANIL 
[86] indicating, in agreement with the coalescence model, that deuterons are on average 
emitted earlier than protons. 

For ultra-relativistic heavy ion collisions, the sensitivity of the correlation 
ratio to the relative time shift (At) (introduced ad hoc) was studied for various two- 
particle systems simulated using the transport codes [85]. The scaling of the effect with 
the space-time asymmetry and with the inverse Bohr radius a was clearly illustrated. It 
was concluded that the 1Z+/1Z- ratio can be sensitive to the shifts in the particle emission 
times of the order of a few fm/c. Motivated by this result, the correlation asymmetry for 
the K + K~ system has been studied in a two-phase thermodynamic evolution model and 
the sensitivity has been demonstrated to the production of the transient strange quark 
matter state even if it decays on strong interaction time scales [87]. The method sensitivity 
to the space-time asymmetries arising also in the usual multiparticle production scenarios 
was demonstrated for AGS and SPS energies using the transport code RQMD [65, 83, 84]. 
At AGS energy, the Au + Au collisions have been simulated and the 7rp correlations have 
been studied in the projectile fragmentation region where proton directed flow is most 
pronounced and where the proton and pion sources are expected to be shifted relative 
to each other both in the longitudinal and in the transverse directions in the reaction 
plane. It was shown [84] that the corresponding TZ+/TZ- ratios are sufficiently sensitive 
to reveal the shifts; they were confirmed in the directional analysis of the experimental 
AGS correlation data [88]. 

At SPS energy, the simulated central Pb + Pb collisions yield practically zero asymme- 
tries for 7r + 7r~ system while, for n ± p systems, the LCMS asymmetries are (Ax) = —6.2 
fm, (Ay) = (Az) = 0, (At) = —0.5 fm/c, (Ax*) = —7.9 fm in the symmetric midra- 
pidity window 8 [83] and, (Ax) = -5.2 fm, (Ay) = 0, (Az) = -6.5 fm, (At) = 2.9 
fm/c, (Ax*) = —8.5, for the NA49 acceptance (shifting the rapidities into the forward 
hemisphere) [65]. Besides, (x) increases with particle p t or u t = Pt/m, starting from zero 
due to kinematic reasons. The asymmetry arises because of a faster increase with u t for 
heavier particle. The non-zero positive value of (x) = (r t x) (x = Pt/pt and r t is the 
transverse radius vector of the emitter) and the hierarchy (x n ) < (xk) < (x p ) is a signal 
of a universal transversal collective flow [65, 66]. To see this, one should simply take into 
account that the thermal transverse velocity p F is smaller for heavier particle and thus 
washes out the positive shift due to the transversal collective flow velocity (3 F to a lesser 
extent. More explicitly, in the non-relativistic approximation, the transverse velocity 
(3 t = (3 F + /3 T ; in the out-side decomposition, f3 t = p t {l,0}, f3 F = p F {cos(f) r ,sm(f) r }, 
(3 T = /Jr{cos0T,sin0<r}- Due to the azimuthal symmetry, the vector of the transversal 
collective flow velocity f3 F is parallel to the transverse radius vector r t = r t {cos0 r , sin0 r } 
and, its magnitude depends only on r t : (3 F = P F (r t ). To calculate (x), one has to average 
over four variables r t , <f) r , Pt and <f) F - At a fixed transverse velocity vector f3 t , only two of 
them (e.g., r t ,(p r or r t , fir) are independent. In particular, (5 T = Pi + Pp — 2p t p F cos(f) r , 

8 (Ay) = due to the azimuthal symmetry and (Az) = in a symmetric mid-rapidity window due to 
the symmetry of the initial system. 
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so the destructive effect of the thermal velocity Pt on the out shift is clearly seen: 

(x) = {r t coa<f> r ) = (r t j . (11) 

The maximal out shift (x) max = (r t ) corresponds to zero thermal velocity. The shift 
vanishes when the width of the contributing interval \(5 t — (5p\ < Pt < Pt + Pf becomes 
negligible compared with the characteristic width of the thermal distribution, e.g., at 
fit or p F — > or, for very light particles; the angle (f> r is then decorrelated from p T 
and so distributed uniformly in the full angular interval (— ir, ir). 9 As a result, in case of a 
locally equilibrated expansion process, one expects a negative asymmetry (Ax) = (x± — x 2 ) 
provided mi < m 2 . Moreover, this asymmetry vanishes in both limiting cases: Pf <C Pt 
and Pf 3> Pt- 

These conclusions agree with the calculations in the longitudinal-boost invariant hy- 
drodynamic model. Thus, assuming a linear non-relativistic transversal flow velocity pro- 
file Pf = Port/r , the local thermal momentum distribution characterized by the kinetic 
freeze-out temperature T and the Gaussian density profile exp(— r 2 /(2rQ)), one confirms 
a faster rise of (x) with p t for heavier particles (see the non-relativistic limit of Eq. (30) 
in Ref. [89]): 

{x) = W+rJmt' (12) 

The maximal magnitude of the asymmetry (x-y — x 2 ) at P u = P 2 t = v t is achieved for an 
optimal value of the flow parameter p = T/(mi t m 2t ) 1 / 2 = T/(7 t 2 mim 2 ) 1 / 2 ; e.g., for np 
pairs at v t = 0.6 (close to a mean LCMS velocity of low-Q 7rp pairs in the NA49 experiment 
at SPS [65]) and T = 120 MeV, the optimal value Po = 0.27. The SPS data on particle 
spectra and interferometry radii in central Pb + Pb collisions at 158 AGeV are consistent 
with the parameters Po 0.35, r ~ 6 fm and T 120 MeV with the uncertainties of 
10 — 20% [20, 35, 90]. The corresponding out asymmetry for np pairs (Ax) = (x^ — x p ) ~ 
—4 fm at Vt = 0.6. As for the longitudinal and time shifts, in the longitudinal-boost 
invariant hydrodynamic model z = rsinh?? and t = t cosh r], where r is the proper freeze- 
out time and r\ is the emitter rapidity. At a given p t , the LCMS ^-distribution of the 
contributing emitters is given by the thermal law exp(— m t coshrj/T). Being symmetric, 
it predicts vanishing longitudinal shift: (z) = (rsmhr]) = 0. To estimate the time shift, 
for m t > T one can write coshr? « 1 +r] 2 /2 and get (t) ~ r(l + |T/m t ). 10 For the central 
Pb-\-Pb collisions at SPS, r ~ 8 fm/c and the relative time shift (At) = (t n — t p ) ~ 3 fm/c. 
This shift is about the same as predicted by RQMD for the asymmetric NA49 rapidity 
acceptance. The magnitude of the relative out shift in pair rest frame (determining the 
observable asymmetry), (Ax*) « —7 fm, is however lower than in RQMD due to ~ 20% 
lower magnitude of (Ax). 

In fact, the NA49 data on 1Z +X /1Z- X ratio for n + p and n~p systems show consis- 
tent mirror symmetric deviations from unity, their size of several percent and the (In- 
dependence being in agreement with RQMD calculations corrected for the resolution and 
purity [66, 70, 72]. Similar pattern of the correlation asymmetries has been reported 

9 Note that, irrespective of the thermal width, the side shift (y) = (rt sin0 r ) = since, due to azimuthal 
symmetry, the angles (f> r and — (f> r contribute with the same weights. 

10 One also recovers the expression for the LCMS interferometry longitudinal radius squared [19]: R 2 Z = 
((z — (z)) 2 ) ~ T 2 T/m t up to a relative correction 0(T/m t ). 
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also for ir ± K ± and systems in experiment STAR at RHIC. They seem to be in 

agreement with the hydro dynamic type calculations with a stronger transverse flow than 
at SPS and a box- like density profile (blast wave), and - somewhat lower than RQMD 
predictions [34, 46]. 

The finite widths of particle rapidity distributions require however a violation of the 
boost invariance. It can be parameterized by a Gaussian dispersion Arj 2 of the LCMS 
^-distribution centered at —Y, where Y is the CMS pair rapidity; e.g., the data on central 
Pb + Pb collisions at 158 AGeV are consistent with Arj = 1.3 [20]. As a result, 

(z) « -tY{\ + Arfmt/T)- 1 (13) 

and (t) acquires a F-dependent contribution |rF 2 (l + Arfm t /T)~ 2 . For the asymmetric 
NA49 rapidity acceptance, the mean irp pair rapidity Y ~ 1.5, (z n — z p ) ~ —2.8 fm 
and the np time shift at Y = is increased by ~ 0.7 fm/c. This is in qualitative 
agreement with the RQMD predictions for the rapidity dependence of the longitudinal 
and time shifts. The magnitude of the F-dependent shifts in the hydrodynamic model 
is however substantially smaller. Besides, the LCMS emission times in RQMD are by a 
factor of 2 — 3 larger and show substantial dependence on the transverse velocity [65]. 
These differences may point to the oversimplified space-time evolution picture in the 
hydrodynamic model. Particularly, the neglect of r r dependence of the proper freeze-out 
time and of the longitudinal acceleration during the evolution may not be justified [20, 89]. 



8 Spin correlations 

The information on the system size and the two-particle interaction can be achieved also 
with the help of spin correlation measurements using as a spin analyzer the asymmetric 
(weak) particle decay [65, 91, 92]. Since this technique requires no construction of the 
uncorrelated reference sample, it can serve as an important consistency check of the 
standard correlation measurements. Particularly, for two A-particles decaying into the 
p7i~ channel characterized by the asymmetry parameter a = 0.642, the distribution of the 
cosine of the relative angle between the directions of the decay protons in the respective 
A rest frames allows one to determine the triplet fraction p t = TZ t /TZ, where lZ t is the 
triplet part of the correlation function (see the footnote in connection with Eq. (8)): 



dN/d cos 9 = l - 



1 + a 2 (j^Pt — cos 9 



(14) 



Both the correlation and spin composition measurements were recently done for two-A 
systems produced in multihadronic Z° decays at LEP [79, 93]. Except for a suppression 
at Q < 2 GeV/c, the triplet fraction p t was found to be consistent with the value 0.75, as 
expected from a statistical spin mixture. Such a suppression, as well as similar suppression 
of the usual correlation function, is expected due to the effects of QS and a repulsive 
potential core, and points to a small correlation radius r < 0.5 fm [80]. 

The spin correlations allow also for a relatively simple test of the quantum-mechanical 
coherence, based on Bell-type inequalities derived from the assumption of the factorizabil- 
ity of the two-particle density matrix, i.e. its reduction to a sum of the direct products 
of one-particle density matrices with the nonnegative coefficients [92]. Clearly, such a 
form of the density matrix corresponds to a classical probabilistic description and cannot 
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account for the coherent quantum-mechanical effects, particularly, for the production of 
two A-particles in a singlet state. Thus the suppression of the triplet AA fraction ob- 
served in multihadronic Z° decays at LEP indicates a violation of one of the Bell-type 
inequalities p t > 1/2. 

9 Conclusions 

Thanks to the effects of quantum statistics and final state interaction, the particle mo- 
mentum and, recently, also spin correlations give unique information on the space-time 
production characteristics and the collective phenomena like multiboson and coherence 
effects and collective flows. Besides the flow signals from single-particle spectra and like- 
meson interferometry, rather direct evidence for a strong transverse flow in heavy ion 
collisions at SPS and RHIC comes from unlike particle correlation asymmetries. Being 
sensitive to relative time delays and collective flows, the correlation asymmetries can be 
especially useful to study the effects of the quark-gluon plasma phase transition. The 
correlations yield also a valuable information on the particle strong interaction hardly 
accessible by other means. 
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